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^ ■ Coherent pion photoproduction, d(j, 7r°)d, is calculated in a coupled channels 

O ' 

^sO ■ approach in momentum space using an elementary production operator that 

Q\ 

includes Born terms, vector meson exchange and M1/E2 A-resonance exci- 
"q ■ tation. The final state interaction is treated in the KMT multiple scattering 

approach that describes ird elastic scattering well. Differential cross sections 
are found to agree well with existing (7, 7r°) data. The polarization observ- 

X 

ables T,,iTu and the tensor analyzing powers are predicted for forthcoming 
experiments. A number of transparent relationships between polarization 
observables and the elementary production operator is derived. Finally, cal- 
culations are performed for coherent eta photoproduction, d{~f, rf)d, in the 

threshold region which agree well with a very recent Mainz experiment. 
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I. INTRODUCTION 



For a long time, coherent pion photoproduction on the deuteron has been studied as 
a source of information on tt° photoproduction off the neutron. In the beginning, the im- 
pulse approximation was the main instrument for this undertaking. However, in contrast 
to charged pion photoproduction where the largest contribution comes from the well-known 
Kroll-Ruderman amplitude, the isoscalar ir° amplitude on the deuteron target is drastically 
reduced due to an accidental cancellation between leading terms. Therefore, the contribu- 
tions from pion rescattering or final state interaction (FSI) become extremely important. 
This dramatic effect from pion rescattering with charge exchange contributions was first 
found by Koch and Woloshyn |l] and then verified by Bosted and Laget || in studies of 
coherent pion photoproduction on the deuteron in the threshold region. 

With increasing photon energy the contribution from the direct term (without rescatter- 
ing) becomes very large due to the dominance of the A-resonance contribution and finally 
dominates the cross section. In the A-resonance region the main mechanism of FSI is elastic 
pion scattering; the contributions from the charge exchange reactions become small. In re- 
cent years, the three-body approach has been employed for the description of FSI. By solving 
the Faddeev equations it is possible to take into account the coupling with the break-up and 
pion absorption channels. In this framework considerable progress has been achieved in the 
description of pion-deuteron elastic scattering (see, for example, the review of Ref. [§]). Very 
recently, Wilhelm and Arenhovel j|,[5| used an approach of NN — NA coupled channels for 
describing deuteron photodisintegration, 7 d — > pn and coherent pion photoproduction, 

-f + d^n° + d. (1) 

In another approach, recently developed for reaction (1) by Garcilazo and de Guerra ||, 
relativistic Feynman diagrams have been evaluated. Using the spectator-on-mass-shell pre- 
scription, differential cross sections and polarization observables have been calculated in the 
energy region from threshold up to 1 GeV. Blaazer |7|, on the other hand, used a Faddeev 
approach to the ttNN system to describe pion rescattering in the final state. 
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This paper is the continuation of our previous work which was devoted to pion scatter- 
ing m| and pion photoproduction [fT0| , |ri}] on very light nuclei. Using a microscopic approach 
based on the KMT multiple scattering approach |12[ in momentum space we have achieved 
a good description of pion scattering on deuteron and 3 He and pion photoproduction on 
3 He. 

In Part I ||, we demonstrated that a microscopic description of elastic ird scattering 
in the framework of the multiple scattering theory is able to describe the differential cross 
section and polarization observables equally well as in the Faddeev approach. Furthermore, 
our multiple scattering approach proved to be very transparent, simplifying a comparison 
of the pion interaction with heavier nuclei. We believe the same rational holds in the case 
of coherent ir° photoproduction. 

Very recently, interest in the physics with eta mesons has grown significantly. This 
relatively new field of theoretical and experimental studies is promising to be very useful 
to obtain a more global picture of the meson-nuclear interaction. On the other hand, these 
studies could give important information on the excitation of nucleon resonances in nuclei 
(like the S , n(1535) and the -Di 3 (1525)). Eta photoproduction on the deuteron is of special 
interest since, due to isospin selection rules a measurement of this process should allow us to 
completely determine the proton and neutron electromagnetic amplitudes and to separate 
the isospin components of the nucleon resonances. 

Unfortunately, despite its long history, the experimental data available for the meson 
photoproduction processes on the deuteron are not of as high a quality as the ones for, i.e., 
elastic pion scattering. However, with the new accelerators at Mainz, Bonn, Saskatoon and, 
in the near future, CEBAF we hope that such data will be available soon. 

The basic ingredients of our formalism based on the multiple scattering approach and 
coupled-channels method are given in Section 2. In Section 3 we present the main expressions 
for the differential cross section and polarization observables. Our results for coherent 7r° 
photoproduction are discussed in Section 4. In Section 5, we will complete our analysis 
with the investigation of coherent eta photoproduction on the deuteron. A summary and 
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conclusions are given in Section 6. 



II. BASIC INGREDIENTS OF THE THEORY 

To obtain the amplitude of pion photoproduction on nuclei one may start, as in pion 
scattering, with the amplitude of the elementary process and the use of multiple scatter- 
ing theory fl2| , [T3|l . An extension of the KMT approach [[TjJ to the case of nuclear pion 
photoproduction leads to the following expression for the T-matrix: 

A A A AAA 

T,(E) = 5>/ + E E * G (E) r/ + E E E ^ G (E) n G (E) r] + ... , (2) 

j=l i^j j=l k^i j^i j=l 

where Tj and tJ describe the pion scattering and pion photoproduction on bound nucleons, 
respectively, and Gq(E) is the Green's function for a free pion- nuclear system. Together 
with the equations for the pion-nuclear scattering T-matrix (see Part I), Eq. (2) can be 
rewritten as a system of integral equations: 

T y (E) = U 7 (E) + T'(E)G (E)U 1 (E) , (3a) 

T'{E) = U\E) + U'(E)G (E)T'(E) , (3b) 

where 

T'(E) = ^T(E), U'(E) = ^j:r J (E),and U y (E) = E r](E) . (4) 

A A 3=1 3=1 

For a consistent derivation of the pion photoproduction amplitude from Eq. (3) we will 
use the same approximations as in pion scattering. First in the spectral decomposition of 
the Green's function only the deuteron ground state is retained, neglecting the contributions 
from the coupling with the break-up channels. As we have seen in Part I, this approximation 
is a reasonable starting point to study pion-nuclear scattering. However, as we will see 
below, near threshold for isoscalar nuclei like the deuteron the ground state contribution is 
very small and the coupling to the break-up channels becomes important. Second we use 
the free elementary t-matrix, ^(E) w r'(uj). The last assumption is called modified impulse 
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approximation where the connection between pion-nucleon (u) and pion-nuclear (E) energies 
is the same as in elastic pion scattering (three-body choice): 

r2 , (q + k) 



u = E + m n + M 



A? 



im^ + M N f+^ +kf 



M 



(5) 
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where m n and M N are the pion and nucleon masses, k and q are the photon and pion 
momenta, respectively. 

Finally, using these two approximations for the nuclear pion photoproduction amplitude 
we get an expression in the distorted wave impulse approximation (DWIA). In momentum 
space it can be presented as 

-r M f Mi{Qi k) — V MfM .{q, K) (27r)2 ^y M{ql) E (q)-E(q')+ie ' W 

where the amplitudes F^ M . and F' M M , are connected with the T 7 (E) and T / (£ T )-matrixes 
in Eq. (3) via relations 



F&lMtf,® = - V " V y < Tr(g),/ | T 7 (£) I i, 7 (a) >, (7a) 



A4(g)A4(fc) 
27 



■^MfM'f ?') = ~ V 2?r < Tr(g), / | T'(E) I /', vr(g') > . (7b) 

In Eqs. (6-7) A = ±1 is the photon polarization, | i >, | / > and | /' > are the 
deuteron ground state which differs only by projection of the deuteron's spin, M iy Mf 
and M'f. The relativistic pion-nuclear and photon-nuclear reduced masses are given by 
M(q) = E 7r (q)E A (q)/E(q) and M(k) = kE A (k)/E(k), where E = E(q) = E(k) is the 
pion-nuclear (or photon- nuclear) total energy in the pion-nuclear (ft A) cm. system. 

Note that the main difference compared to the standard DWIA approach (which is not 
appropriate for the deuteron) is the factor (A — I) /A in the T'-matrix. This factor avoids 
double counting of pion rescattering on one and the same nucleon. Such effects are already 
included in the elementary amplitude. Another difference is that Eq. (6) includes not only 
coherent (non spin-flip) pion rescattering but also incoherent spin-flip transitions which 
change the projection of the spin of the deuteron. 



In the framework of the plane wave impulse approximation (PWIA) the pion photopro- 
duction amplitude is equal to the first term of Eq. (6) which is expressed in terms of the 
elementary pion-nucleon photoproduction t 7 (o;)-matrix or amplitude 



m - JM(q)M(k) * 

J Z7T ■ , 

J = l 

< «f,p" | *» | p, k\ >= - ^ / (A) (^, C) <*(p' + 9 - P - *) , (8b) 

^(q,p')fM{k,p) 

where relativistic pion-nuclear and photon-nuclear reduced masses are fi{q,p') = 
E n (q)EN{p')/uj and fi(k,p) = kE^ip)/^, p and p' are the nucleon momenta in the initial 
and final states (in it A cm. system), respectively. 

The relativistically invariant amplitude f( x \ in the pion-nucleon (ttN) cm. system has 
the standard CGLN form 

f w {uj, 6*) = iF t a -e x + F 2 a^qa-[kx e x ] + iF 3 a ■k'q-e x + iF 4 a ■ "q"q- e x , (9) 

where q and k are the unit vectors for the pion and photon momenta in the nN cm. system, 
8* is the pion angle in the same system. All these variables can be expressed in an arbitrary 
frame via the standard Lorentz transformation (see Ref. [|TT|1 ). 

In nuclear applications, it is convenient to divide the amplitude in Eq. (9) into non 
spin-flip (/ 2 ) and spin-flip (fx, f 3 , / 4 ) amplitudes 

fW = if x 3 ■ e x + f 2 [k x q] ■ e x + if 3 a -k q-e x + ifaa -qq-e x . (10) 

These are connected to the standard CGLN amplitudes via 

f 1 = F 1 -F 2 cos6l, f 2 = F 2 , f 3 = F 3 + F 2 , / 4 = F 4 . (11) 

In terms of the lowest s- and p-wave multipoles, the amplitudes fi can be expressed as 

fx = E 0+ + cos 0^(3 £ 1+ + M 1+ - M x _) , (12a) 

/ 2 = 2M 1+ + M 1 _, (12b) 
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f 3 = 3E 1+ -M 1+ + M 1 _, (12c) 



U = 0. (12d) 



In our calculation, we use the unitary version of the Blomqvist-Laget amplitude [ ]IoT , |Tt3| for 



the elementary amplitude which contains contributions from Born and cu-exchange terms. 
The dominant A-resonance contribution is described by the real and imaginary parts of 
the resonant M^ + {uj) and E^ + (uj) multipoles that are added to the background terms in a 
unitarized way. In the framework of this model the elementary process is well described for 
photon energies up to i? 7 =500 MeV. 



The explicit expression for the full elementary amplitude is given in Ref. ||11|| . We only 
remind that in accordance with the Blomqvist-Laget approach it is obtained after a non- 
relativistic reduction (up to order (p/Mjv) 2 ) of the relativistic and gauge invariant amplitude. 
One of the advantages of this model is that the final expression for the elementary amplitude 
is well appropriate for nuclear calculations with non-relativistic wave functions. However, 
due to the reduction procedure the gauge invariance is lost. The neglected terms that could 
restore the gauge invariance are of the order of (p/M^) 3 (consistent with our approximation) 
and their influence is small. 

Finally, we mention the so-called factorization approximation in the treatment of nucleon 
Fermi motion. We have developed the formalism for the exact treatment of the nucleon 
momentum dependence in the elementary amplitude using a multidimensional integration 
in Eq. (7). However, similar to our previous studies [0,0 the substitution 



P^Peff. = -\- ^/(k ~q) = ~-\0 (13) 

works very well quantitatively and provides an excellent approximation. This result is based 
on the fact that in case of a Gauss wave function, which roughly reproduces the dominant 
S-wave part of the deuteron ground state, the replacement (13) gives an exact treatment for 
the linear p/2M terms in the elementary amplitude. 
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III. POLARIZATION OBSERVABLES 



Polarization observables have the promise of opening a new field in the electromagnetic 
production of pions from protons and nuclear targets. Since many of these observables 
contain interference terms, small but important amplitudes can be investigated in a unique 
way. 

To define our polarization observables we choose the system of coordinates corresponding 
to the Madison Convention ||19|| : a right-handed coordinate system where the positive z-axis 
is along the photon beam direction and the y-axis is along the vector [k x q\. 

In this paper, we focus only on single polarization observables that appear in ir° photo- 
production with a polarized beam and an unpolarized target, or an unpolarized beam and 
a polarized target. In both cases five polarization observables can be measured: 
The photon asymmetry 

da/dQ^ + da/dfl(W) ' ^ ' 

where _L (||) refers to a photon linearly polarized perpendicular (parallel) to the reaction 
plane and four analyzing powers (the vector, iTn, and three tensor ones, T2o,T2i,T 2 2) 

Tr(Fr kK Fl) 

Tr(FFt) ' 1 J 

In Eq. (15), the spherical tensor operator t^ k is defined as 

< J M' I r kK I JM >= Jk(-l) J+M ' ( ^ \ t k ) (16) 



with a Wigner 3j symbol and J = y/2J + 1. 



From the partial wave decomposition of the electromagnetic vector potential 13 it 



follows that the pion photoproduction matrix elements, F^] M ., can be divided into magnetic 
and electric components which are both independent of the photon polarization A = ±1: 

p (A) _ p (M) , p (E) 

P M s Mi — " MjMi + A ^M f Mi- 



s 



Similar to the formalism in pion scattering, one can present the photoproduction matrix in 
the following form 

\ 



( jpW p(A) 771(A) \ / 
r ++ r 0+ r -+ 



? (A) 
M s M l 



r (A) P (A) 
+0 -^00 



i(A) 
-0 



V 



.(A) Fq (A) ^(A) 



V 



Ae Be Ce x 



+ A 



/ 



V 



V E X? s 
-C E Be -A e 



(18) 



Bm Cm 
Cm —Bm Am 

where the signs +,0,— correspond to the deuteron spin projections = +1,0,-1, 

respectively. Note that the structure of the magnetic part of matrix (18) is similar to the 
Robson matrix which was discussed in the case of elastic pion scattering ||. The main 
difference appears in the electric part, whose contribution is much smaller. Therefore, the 
behavior of the differential cross section and the single polarization observables is mainly 
determined by the five magnetic amplitudes. We point out that the Ce amplitude vanishes 
in PWIA; it has contributions due only to pion rescattering. 

Calculating the expectation values of the operator r^ K for k = k = 0, the differential 
cross section and photon asymmetry can be written as 



da doM da E 

In = ~dh~ + In 



(19a) 



da M _ 2q 2 

In ~ Am 1 



Cm | 2 + I T^m | 2 +— I &m | 2 ) 



(19b) 



da E 2q 
lh~ = 3k { 



A E \ + \ B E \ + \ C E \ + \ V L 



daM I dfl — daE / dn 



(19c) 



(19d) 



daM/dVt + das/dVL 

Similarly, the vector and tensor analyzing powers can be divided into magnetic and 
electric parts 



rp _ rp(M) T (E) 
J-kn — J-kK i 1 kn 



(20) 



The expression for T^f' in terms of the Am,---,£m amplitudes are the same as in pion 
scattering: 



iTlP = ^-Im[V M {A M - C M ) - B* M £ M ) }/a(9 n ) , (21a) 

T 2 ( M) = ^(| Am I 2 + I B M | 2 + | C M | 2 -2 | V M | 2 - | £ M DM^) , (21b) 

^2 ( i M) = ~^Jl Re \ V *m (Am ~ C M ) + B* M E M ) ]/a(9 n ) , (21c) 

T 2 ( 2 M) = -L[2Re(A M C M )- | #m | 2 ]M^), (21d) 

. „ . k da . 

«.<«.) = (21.) 

For the electric part of the analyzing powers we obtain 

^if ) = ^Im[V* E (A E + C B ) ]/a(^) , (22a) 

Tio ] = ^ (I ^ I 2 + I B B | 2 + | C E I 2 -2 | D E | 2 )/a(^) , (22b) 

T£ ] = -^Re[V E (A E + C E ) }/a(9„) , (22c) 

T 2 ( 2 E) = -J-[2Re(A* E C E )- I 8* | 2 ]/a(^) . (22d) 



Finally, we present expressions for the A, B, etc. amplitudes obtained: 1) in the plane 
wave impulse approximation, 2) by taking into account only contributions from the S- 
component and its interference with the D-component of the deuteron wave function. Then 
the Robson amplitudes can be expressed solely through the elementary amplitudes fi, ft 
from Eq. (10) and the radial integrals R^ (introduced in Part I): 

RfhQ)= j drUvWM^Ufr), (23) 

where Ui is the deuteron radial wave function Jl{z) is the spherical Bessel function, and 
{/, V} = 0, 2. This yields 
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Am = V2f 2 sin 91 



d(0) 
-"-00 



2V2 



(1 + 3 cos 29) R t 



(2) 
02 



w A , 



(24a) 



B M = [fi Rs + ^=/ 2 sin 0; sin 29 4? j , 



(24b) 



C M = -3/2 sin 0; (1 - cos 29) i$ W A , 



= - (/1 i?5 - ^=/ 2 sin 61 sin 29 , 



(24c) 



(24d) 



£ M = V2f 2 sm9l 



4° ) + - ? =(l + 3cos29)4 2 2 ) 



W A , 



(24e) 



^ £ -v^sin^ 



(/a + U cos 0;)i2 s + -j=f Q cos 9 i? ( 



(2) 
L 02 



(24f) 



(/1 + / 4 sin 2 e;)R s - -t=/q sin 0; sin 9 i2, 



V2 J 



(2) 
02 



w A , 



where we have used the following notations: 



(24g) 



Rs(Q) = R§}(Q) + ^R$(Q), 



(25a) 



fq = + /a cos 0; + / 4 ) - |(/ s + / 4 cos 01) , 



(25b) 



cos 9 = —(k — g cos 6^) . 



(25c) 



— * — * 

In Eqs. (24-25) #* is the pion angle in tcN cm. system, Q = k — q is the nuclear momentum 
transfer. The kinematical factor W A can easily be obtained from Eq. (8). 

As we will see below, the expressions (24) are very useful for a qualitative as well as 
a quantitative understanding of the polarization observables. The complete expressions 
(including complete contributions from the deuteron D-state component) are given in the 
Appendix. 
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IV. RESULTS AND DISCUSSION 



A. Differential cross section 

We begin our discussion with the results for the differential cross section. Using our 
formalism described above we will show that it reproduces the well-known results for the 
coherent tt° photoproduction at threshold and in the A-resonance region. 

In general, pion photoproduction near threshold is dominated by the spin-flip amplitude 
with the leading multipole E 0+ . However, neutral pion photoproduction is comparable to 
the situation with the scattering length in elastic pion-deuteron scattering (see Part I): the 
isoscalar part of the s-wave multipole is one order of magnitude smaller than the isovector 
one. Therefore, the PWIA approach is not valid at threshold and contributions from pion 
rescattering with charge exchange become very important |l],j2|. For this two-step process 
the corresponding contribution to the total amplitude can be estimated by using the first 
iteration of Eq. (3) 

41. = — W A < d | e* -^ e — £ f„ N (i) / W (j) | d > , (26) 

where Ai^d and p, n N are the reduced masses for the nd and 7iN systems, respectively, | d > 
is the deuteron ground state and go is the on shell momentum of the pion in the intermediate 
state. 

In the limit of q — > 0, the elementary amplitudes for ttN scattering, f n ^, and pion pho- 
toproduction, /W, become constant and the total amplitude (with direct and rescattering 
contributions) can be reduced to the expression |^TJ 



F^^W A R { S(Q) 



i L-""l' A'~"'M Mird j 1 \ o n + , + n k p ^ 

(£ + + E 0+ > H \ - / a 7rN l^o+ - A 

H'ttN \ ' 



(27) 



with the pion-nucleon scattering length 0%$ = —a%N = V^( a i ~~ a 3)/3 = 0.122/m^ 



and < 1/r >= 0.65/m w . In Eq. (27), the Eq+ n multipoles for ir° photoproduction are 



much smaller than the Kroll-Ruderman multipole E£ + N , thus the rescattering terms are 
comparatively enhanced. 
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In Fig. 1 we show our calculations for the total and differential cross section in the 
threshold region. In this region the contribution of elastic rescattering (without charge 
exchange) is very small. This is due to the small size of the isoscalar part of the nN scattering 
amplitude. On the other hand, due to its large isovector part and also due to the large E^ N 
multipole, the role of the charge exchange process becomes very important. In the upper 
figure the reduced cross section is compared with a recent re-analysis of an earlier Saclay 
measurement. The full calculations are about 15% larger than the experimental data. The 
same results have recently been obtained by Garcilazo and Guerra || using a Feynman 
diagram method with a pure pseudovector ttNN coupling. Calculations by Bosted and 
Laget are about 15% lower than ours and agree better with the experiment. 

The authors of Ref. argue that the main reason of their disagreement with the ex- 
periment could be due to the contributions from third- and higher-order scattering terms. 
However, our calculations (compare dashed and solid curves in Fig. 1) show that such contri- 
butions usually increase the cross section because of the attractive nature of the pion- nuclear 
interaction in the threshold region. 

In the lower part of Fig. 1 we show calculations for the differential cross section at 
AE = Ej — E 1 ^' = 8 MeV. Here we see that the interference between direct and rescattering 
terms is strongly angle dependent, leading to an increase in the differential cross section by 
about a factor of two at backward angles. 

Recent developments in ir° photoproduction on the proton p3|j2^] that indicate large 



deviations from the old Low Energy Theorem (LET) predictions suggest that such effects 
could also be visible in coherent photoproduction on the deuteron. In the direct term 
only the isospin (+) amplitude appears which is proportional to the sum of the proton 
and neutron amplitudes, Eq+ = (E 0+ (ir p) + E 0+ (7r°n))/2. While PV Born terms predict 
this amplitude to be constant in the threshold region with a value close to the old LET, 
E$ = -1.0 • KT 3 /m^, the new experiments on the proton show an energy dependent 
function with a cusp effect at the charged pion threshold. Using isospin symmetry, with 
this new experimental data and the older values for (7, ir^) |25| we find values for from 
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-1.5 to +1.5 in units of 10~ 3 /m n . Predictions from chiral perturbation theory give positive 
values for this amplitude that vary between the neutral and charged thresholds in the range of 
(0.5 — 1.2)-10 _3 /m 7r ||26|| . In Fig. 2, we show the influence of this isospin (+) amplitude on the 
total and differential cross section of coherent pion photoproduction on the deuteron. Even 
though the direct term is not dominant in this reaction, a significant sensitivity can be seen 
because of the interference with the pion rescattering channel. Comparing our calculations 
with the total cross section data in Fig. 2a leads to a preferred value of —0.5 • 10 _3 /m 7r . As 
demonstrated in Fig. 2b, a determination of Eq+ should be much easier using angular 
distributions, where one can easily distinguish different values in the forward-backward 
angular asymmetry. 

With increasing photon energy the contributions from p-wave multipoles become 
stronger, due mostly to the excitation of the A-resonance. Fig. 3 illustrates the role of 
the individual parts of the elementary photoproduction amplitude: the Born+cu-exchange 
terms (dotted curves), the A-resonance contribution (dashed curves) and the total amplitude 
(solid curves). The Born and tu-exchange diagrams are essential at forward and backward 
angles. In this region only the f\ amplitude remains which contains the nonresonant E 0+ 
and M]_ multipoles as well as the resonant and amplitudes. Around 8^=90°, the 
A-resonance contribution dominates due to the magnitude of the coherent contribution from 
the non spin-flip amplitude f^- 

In comparison to the available data the overall agreement is satisfactory. An exception 
may be the large disagreement in the forward direction with the data from Ref. |30| with 
rather large error bars. The best description of the experimental measurements is achieved at 
9^ = 90°, where the A contribution dominates. Together with our earlier results obtained for 
3 if e and p-shell nuclei, this gives us confidence in the treatment of A propagation in nuclei. 
In order to draw further conclusions more precise experiments are needed, in particular at 
extreme angles where the A contribution is weaker. 

Fig. 4 shows the differential cross sections at E 1 = 260 — 400 MeV. Here we illustrate 
the role of the spin-flip contribution (dash-dotted curves). It only contributes significantly 
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at very backward and forward angles. Around 8 n = 90° the main contribution comes from 
the non spin-flip transition associated with the elementary amplitude fi which is dominated 
by the A-resonance. 

The role of the rescattering mechanism can be seen by comparing the simple PWIA 
(dashed curves) with our full calculations (solid curves). Its contribution is biggest in the 
A-resonance region and at backward angles. Among the different rescattering mechanisms 
elastic scattering (without charge exchange) dominates in this energy region. This is in 
contrast to the threshold region where , as we have seen above, charge exchange rescattering 
gives the main contribution. Above the A-resonance region, the rescattering contribution 
becomes smaller with increasing photon energy. 

The role of the 3 Di configuration in the deuteron wave function is also illustrated in 
Fig. 4 by turning off the nuclear matrix elements with L = 2. The D-w&ve contribution 
becomes manifest at large angles, 9 V > 90°, for energies E 1 > 300 MeV, which corresponds 
to momentum transfers of Q 2 > 3.1 fm~ 2 . 

We conclude our analysis of the differential cross section by comparing with previous 
calculations. Our results presented here are very similar to the ones obtained by Bosted and 
Laget . The main difference between our and their approach is the treatment of the pion 
rescattering contributions. In our formalism the coupling with the break-up channels was 
not included in contrast to Ref. ||. Apparently, the corresponding contribution is small in 
the kinematical range considered here. 

A comparison with the Garcilazo and Guerra || calculations shows that in the forward 
direction their results are generally about 50% larger, except at £^=340 MeV. This could 
indicate a difference in the spin-flip part of the elementary amplitude which dominates 
in the forward direction. Below we show that this leads to a dramatic difference in the 
vector analyzing power. The work by Blaazer includes only a comparison with the data 
at £"^,=300 MeV. Their results are about 20% above the measured cross section. Within 
their Faddeev formalism they find only small pion rescattering effects. Compared to the 
calculations by Wilhelm and Arenhovel || our results are lower and agree better with the 
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data, except for 6 W = 6°. The main reason for this difference is due to pion rescattering in 
our approach which reduces the cross sections considerably in the region of the A-resonance 
(see Fig. 4). 



B. Polarization observables. 

Photon asymmetry. Fig. 5 shows the angular distribution for the photon asymmetry 
calculated at £ 7 = 200 - 400 MeV in PWIA (dashed curves) and DWIA (solid curves). 
The dotted curves are DWIA calculations without the D-state configuration in the deuteron 
wave function. Clearly, all calculations are very similar; thus, the influence of the deuteron 
.D-state and pion rescattering is small. Moreover, our analysis showed that the contribution 
of the / = 2 multipoles in the elementary production amplitude is also negligible, e.g., we 
can set f± = 0. Thus, using Eq. (24) the photon asymmetry can be expressed in a rather 
simple form using only the elementary amplitudes /i, ji and f^- 

^(deuteron) « \ sin* /t |, j^f^' '+|| A ■ < 2 «) 

This expression is close to the expression for the elementary photon asymmetry for the 
process on the nucleon: 

Unudeon) » \ sin* 0, } - ? -0^1^- j - - . (29) 

Therefore, the angular distribution of the photon asymmetry for the production process on 
hydrogen and deuteron looks very similar. 

Tensor analyzing powers. Fig. 6 presents our results for the T 20 observable obtained in 
the same way as the photon asymmetry. However, in contrast to the photon asymmetry, 
T 2 o is very sensitive to the deuteron D-state component (compare solid and dotted curves). 
The influence of pion rescattering still remains small. 

In the region of 90° < 9 W < 150°, where the contribution of the D-state is more important 
the value for T 20 can be estimated from Eq. (24) as 
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r » ( ~ 130 °' » I a"' H\eJr » - (1 + 3ros29) i + 3^29 ■ (30) 

where angle is defined in Eq. (25c) and x = (Q)/Rqo (Q) ■ This ratio describes the 
role of the deuteron D-state as a function of the momentum transfer Q. It is clear that with 
increasing Q the ratio x also increases and | T 20 | reaches its maximum value. 

The contribution of the deuteron S-state dominates at forward and backward angles. In 
this region Em = Dm = —Be = —T>e and other amplitudes are zero. Therefore, independent 
of the photon energy we have 

T 20 (0°)^T 20 (180 )^-^. (31) 

As can be seen in Fig. 6, this relation is generally fulfilled and only small deviations appear 
due to pion rescattering and Fermi motion. 

Thus, we find that the maximum value of the tensor analyzing power T 20 is only sensitive 
to the .D-state component of the deuteron wave function and not to the details of the 
elementary amplitudes or to pion rescattering. 

The behavior of the other tensor polarization observables T 21 and T 22 is similar (see Fig. 
7 and Fig. 8). They both have a maximum negative value around 9 n = 90°. As in the case 
of T 20 this value is determined by the deuteron D-state. Using Eq. (24), this dependence 
can be expressed approximately as 

r 21 (~90°)^-V6sin 2 e i + 3 J s . n22e , (32a) 

r,(~90V-/|(l-cos2e) 1+3 J s . n . 22e . (32b) 

From this equation it follows that at the maximum T 2 i m 2T 22 . In the forward and backward 
direction these observables vanish. At low energies the S'-state component becomes visible 
and the D-state contribution is negligible. In this region, T 21 can be estimated by the 
following expression 

r<f > (El ~ 200MeVO „ S in ^ - ? +|si ,^ l |2 +| , % p) ■ (33, 
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Vector analyzing power. In our study of polarization observables in elastic pion-deuteron 
scattering (see Part I) we found that the vector analyzing power iT n is very sensitive to 
the details of the theory. Especially pion rescattering was very important, around the A- 
resonance region it changed the sign of iTn. 

Our calculations, shown in Fig. 9, indicate that the influence of pion rescattering is not 
important for iTn in the case of 7r° photoproduction (compare solid and dashed curves). 
The contribution of the deuteron D-state is also small. Therefore, in this reaction the 
vector analyzing power is neither sensitive to details of the rescattering mechanism nor to 
nuclear structure. Therefore, it can be expressed entirely through the elementary amplitudes 
fi, f±. If we again neglect the small contribution from the / 4 amplitude iT u has a simple 
form, 

,^f»(^) W -^, |/i|a+ ff^- a ^| />|r (34) 
Furthermore, due to the small nonresonant E 0+ multipole, f\ reduces to 

h « cos0 w (SE 1+ + M 1+ - Mt_) . (35) 

Therefore, iTu has a sin 9 n cos 9 n angular dependence and changes the sign around 9 n = 
90°. Note that this result is in strong disagreement with the predictions of Garcilazo and 
Guerra || . Their vector analyzing power is positive in the A- resonance region at all angles. 
This can occur only through differences in the elementary amplitude. At the same time 
their predictions for the photon asymmetry and tensor analyzing powers are very similar to 
ours. 

Another interesting consequence which comes from Eqs. (34) and (35) is related to the 
contribution of the Ei + multipole. It is of a particular interest because of its resonant part 
due to the E2 transition in the jNA vertex. The size of this quadrupole component is 
very sensitive to the tensor force in the quark-quark interaction. Thus, this transition could 
give us a measure for the deformation of the delta and, subsequently, of the nucleon in the 
framework of quark models. 

18 



Fig. 10 illustrates the sensitivity of the vector analyzing power to the contribution of 
the E2 transition in 7 A" A vertex at 9 n = 30°. As follows from Fig. 9, the influence of 
other theoretical uncertainties is minimal in this region. However, the presence of the 
multipole dramatically changes the value of iT\\. 

Note that a similar situation also occurs in 7r° photoproduction on the nucleon. Again, if 
we neglect the small contribution from the f± amplitude we obtain for the target asymmetry 
on the proton 



V. ETA PHOTOPRODUCTION 

We complete our analysis with an investigation of coherent eta photoproduction on the 
deuteron. The formalism for this process can be developed in a straightforward way by 
the same coupled channels method which has been applied to pion photoproduction. In 
momentum space the nuclear photoproduction amplitude can be written as 

F m (ff, k) - Vrn (q, k) {2n)2 Mi{q)) £M _ £M>) + . £ , (37) 

where k is the photon, and q is the eta or pion momentum. The total energy in the 7/-nucleus 
and 7r-nucleus channels is denoted by £i(q) = Ei(q) + E^q), the reduced mass is given by 
MM) = Ei{q)E A {q)l£i{q) and a = {A - I) /A. 

V m is expressed in terms of the free eta-nucleon photoproduction t-matrix 



V$(q,k) = - V vK l J < V (q)J\ Strait) I hl(kX) > ■ (38) 

2tt j=1 

As an elementary t vl amplitude we will use the extended dynamical model of Bennhold 
and Tanabe 0,^] where the resonance sector is constrained by solving the coupled channels 
problem for the irN — > ttN, ttN — > rjN, ttN — > ttttN and — > ttN reactions, using available 
data as input. The background is formed by the contributions from s- and li-channel Born 
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terms and by p- and u- exchange in the t-channel. Following the analysis of Ref. [3^] we use 
pseudoscalar Born terms with an rjNN coupling constant of g^/47r = 0.4. 

At present, little is known about the nature of the eta-nucleus interaction which enters 
in the amplitude F V i(q,q'). For the elementary ttN — ► i]N process the experimental data 
is much less complete and accurate in contrast to 7riV-scattering data. There are only few 
theoretical studies of this reaction [ 3l|j33j34 | based on the coupled channel isobar model 



for the nN, rjN and 7T7riV systems. More detailed information about the structure of the 
eta-nucleus interaction can be found in our previous work |§. We only note that, just as 
in pion photoproduction, the eta scattering amplitude F vi is constructed as a solution of 
the Lippmann-Schwinger equation using the KMT version of multiple scattering theory. 
At present our calculations have been carried out in DWIA without coupling to the pion 
channel because our elementary pion photoproduction amplitude is not well appropriate for 
the photon energies E 1 > 500 MeV where higher nucleon resonances become important. 



From the other hand as it was estimated in Ref. |57| contributions coming from the pion 
channel are much smaller than eta rescattering terms. 

In Fig. 11, we show the differential cross section for coherent eta photoproduction on 
the deuteron at E 1 = 675MeV (48 MeV above threshold). The reaction is completely 
dominated by the resonance contribution, the non-resonant Born terms and vector meson 
contributions are negligible. The eta rescattering contribution in the final state is also small 
in this energy region. 

Fig. 12 presents the energy dependence of the differential cross section at 9 = 90°. 



Old experimental data |35j are in dramatic disagreement with numerous theoretical pre- 
dictions [|36| -|38"11 ■ This is related to the small isoscalar (7,77) amplitude E?o+ with R = 
E$/E$ = 0.22. As can be seen in Fig. 12, only an unrealistically large isoscalar ampli- 
tude of Eq_1/Eq+ = 0.7 can explain the data. This puzzle was a subject of investigations by 
many authors over the last twenty years. Recently, a new experiment at Mainz by Krusche 
et al. [39| has obtained an upper limit for coherent rj photoproduction which is in agree- 



ment with our calculations using a realistic ratio of R ~ 0.2. This is also confirmed by a 
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preliminary analysis of an experiment at Bonn 



PI 



VI. CONCLUSION 

In this paper we have presented a calculation for coherent pion and eta photoproduction 
from the deuteron using an approach that is consistent with the description of pion elastic 
scattering (Part I) || . The deuteron wave function used in our computation is obtained from 
the Paris potential. The final state interaction of the emitted mesons is treated in multiple 
scattering as described in Part I. All calculations are performed in momentum space treating 
dependence of the elementary amplitude from mesons and nucleon momenta exactly. 

Our calculations agree well with existing data of differential cross sections and can also 
reproduce previous calculations found in the literature. In the case of tt° photoproduction in 
the threshold region we have found significant sensitivity of the differential cross section to 
the E$ isospin amplitude of the nucleon that should allow a determination of the E 0+ (7r°n) 
amplitude pifl . For coherent rj photoproduction new experimental data from Mainz is in 



agreement with our results obtained using a realistic ratio of E$/E$ = 0.22. 

We have applied our approach to polarization observables that can be measured with the 
new electron scattering facilities where polarized targets and beams are becoming available. 
A number of simple relationships between polarization observables and the basic production 
amplitude was derived for the kinematical regions where nuclear structure and pion rescat- 
tering are not important. We found a strong sensitivity to the deuteron D-state component 
in the tensor analyzing powers. The vector analyzing power (iTn) or target polarization 
shows an enhanced sensitivity to the resonant E^ + multipole. This very important quantity 
may give insight into the medium modification of baryons in the presence of other nucleons. 
A study of this amplitude is possible by looking at different nuclei. In a previous paper on 
3 He it was found that this multipole can be rather easily observed with a polarized photon 
beam. 

Up to now there are only few good data available for this rather elementary reaction 
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and polarization observables are missing in the low energy and A region. With the new 
accelerators at Mainz, Bonn, Saskatoon and, very soon, CEBAF and the development of 
very powerful spectrometers and detectors we hope that soon precise data will become 
available, thus improving our knowledge of the production and interaction of mesons with 
very light nuclei. 
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APPENDIX: 

In this appendix we give PWIA expressions for the Robson amplitudes A, ...£ with the 
deuteron S- and -D-states taken into account. They are expressed with nuclear matrix 
elements Mslj{Q) defined in Ref. || 



M 000 (Q) = 2 ^ RSS(Q) + B$t(Q) 



>(0), 



(Al) 



M 022 (Q) =4 




R { S{Q)-^R < ikQ) 



(A2) 



M W1 (Q) 




(A3) 



M 121 (Q) 



1 47T 



(A4) 



For the magnetic part of the A, ...,£ amplitudes in the PWIA approach we obtain 



A 



M 



'—k sin 61 
6 J - 



M 000 (Q) - -^=(1 + 3 cos 20) M 022 (Q) 



W A , 



(A5) 
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B 



1 4vr 



M 




h M S (Q) + -f 2 sin e* n sin 29 M 022 (Q) 



W A , 



(A6) 



Q I A 

Cm = -g V y/2 sin^ (1 - cos26) M 022 (Q) W A , 



(A7) 



V 



M 



1 /4vr 

2 V ~6~ 



fiM s (Q) - -h sin 0; sin 20 M 022 (Q) 



^A, 



(A8) 



Air 

£m = W— / 2 sin6»* 



Mooo(Q) + — =(1 + 3 cos 26) M 022 (Q) 



2^ 



W A 



For the electric part we obtain 

1 /47T 



2V 3 Sin ^ 



(/a + / 4 cos 0;)M s (Q) - -^/o cos 9M m (Q) 



(A9) 



W^a, (A10) 



1 4tt 




2V 6 



(/1 + A sin^ e;)M s (Q) + -j=f Q sin 0; sin 6 M 121 (Q) 



(All) 



In the PWIA approach the electric amplitude Ce = which contributes only in the presence 
of pion rescattering. In Eqs. (A5-A11) we used the notation 

1 



M S (Q) = M m (Q) - ^=M 121 (Q) 



(A12) 
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FIGURES 

FIG. 1. (a) Reduced total cross section as a function of the photon lab. energy near threshold. 
The dotted curves are the PWIA results, the dashed and solid curves are the calculations including 
pion charge exchange without and with elastic pion rescattering, respectively. The data are from 
Ref. p2[ . (b) Reduced differential cross section for the photon excess energy in the lab. frame 

? 7 - 



above threshold, AE = Ej — E* hr ' = 8 MeV. The notation of the curves is the same as in (a). 



FIG. 2. (a) Sensitivity of the reduced total cross section on the isospin (+) amplitude in the 
direct term. From top to bottom the Efu) s-wave amplitude has been taken as —1, —0.5, 0, +0.5, +1 
in units of 10 _3 /m vr . All curves have been calculated with full rescattering including the SCE 
contribution. The data are from Ref. |]22[. (b) Reduced differential cross section for the photon 



excess energy in the lab. frame above threshold, AE = E^ — E^ hr ' = 8 MeV. The notation of the 
curves is the same as in (a). 

FIG. 3. Energy dependence of the differential cross section at fixed angles 6 W = 6°, 90°, 120° 
and 175°. The dotted and dashed curves are the calculations only with Born+w-exchange and A 
contributions, respectively. Solid curves are the results of the full calculations. Experimental data 



are from Refs. |27| (o), f2§ (•) and gg] (A). 



FIG. 4. Angular distribution at fixed photon lab. energy £ 7 = 260, 300, 340 and 400 MeV. The 
dashed and solid curves are PWIA and full calculations. The dash-dotted and dotted curves are 
the calculations without non spin-flip and without deuteron L>-state, respectively. Experimental 
data are from Refs. |7| (o), Refs. p| (•) and @ (A). 



FIG. 5. Photon asymmetry £ at £^=200-400 MeV. The solid and dashed curves are full and 
PWIA calculations, respectively. The dotted curves are obtained without deuteron D-state. 

FIG. 6. Tensor analyzing power T20. The notations of the curves are as in Fig. 5. 

FIG. 7. Tensor analyzing power T^\- The notations of the curves are as in Fig. 5. 
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FIG. 8. Tensor analyzing power T22. The notations of the curves are as in Fig. 5. 



FIG. 9. Vector analyzing power T\\. The notations of the curves are as in Fig. 5. 

FIG. 10. Sensitivity of iT\\ to the multipole. The solid and dashed curves are the results 
obtained with and without this multipole, respectively. 

FIG. 11. Differential cross section for coherent eta photoproduction on the deuteron at 675 
MeV photon lab. energy. The dotted curves show the individual contributions for the resonance 
excitation and the non-resonant background from vector mesons and Born terms in PWIA. The 
dashed and full lines give the full calculations in PWIA and DWIA, respectively. 

FIG. 12. Energy dependence of the differential cross section at fixed angles 6 V = 90° for the 
d(pf,rj)d reaction calculated with different ratios R = Eq+/Eq+ = 0.22, 0.5, and 0.7, where the 



first value is the prediction of the Bennhold Tanabe model |31|. The dashed and solid curves are 



the PWIA and DWIA calculations, respectively. Experimental data are from Ref. [39 (o) and [35] 
(•)• 
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